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ABSTRACT 


'l  l>c  iiatui;il  convection  heat  tiansfcr  cliaiactciistic*:  of  a  3x3  array  of  vertically  ori- 
cntcil  heated  pi oti nsion*:,  iininersed  in  a  dielectric  lir|nid,  were  investigated.  Aluininuin 
blocks,  2dmin  \  Sinin  x  6niin,  were  used  to  sinuilatc  20  pin  dual  in-line  packages.  Sur¬ 
face  temperature  measurements  of  the  components  were  made  by  imbedding  copper- 
constantan  thermocouples  below  the  surface  of  each  component  face.  A  constant  heat 
Mux  was  piovidcd  to  each  component  using  an  Inconel  foil  heating  element.  Power 
supplied  to  each  component  vaiied  from  0. 1 15W  to  2.90W.  I  he  aluminum  blocks  were 
mounted  on  a  plexiglass  substiate  to  form  a  3x3  array  of  simulated  electronic  compo¬ 
nents.  I  he  ciicnit  boaid  containing  the  components  was  placed  in  a  rectangular, 
plexiglass  enclosure  with  inner  dimensions:  L=  203.2mm  11  =  1.52.()min  W=  82.6mjn.  and 
a  wall  thickness  of  2.3.4111111.  1  he  upper  boundaiy  was  maintained  at  lO^'C,  while  all 

other  exteiior  surfaces  weie  insulated,  I  he  chamber  width,  measured  from  the  surface 
of  the  circuit  boaid  to  the  opposite,  inner  wall  of  the  enclosure,  was  varied  from  42mm 
to  7mni  In  inserting  plexiglass  spacers  into  the  enclosure.  I  wo  dielectric  fuiuids,  rC-75 
and  rc%43,  were  used  as  woiking  fluids. 

Non-dimensional  data  from  this  study  was  combined  with  the  data  obtained  by 
A\tar  (1991)  for  a  hoii/oiital  component  oiicnlation,  to  devclope  an  emperical  corre¬ 
lation  which  predicts  the  Nusscll  number  as  a  function  of  Rayleigh  number.  Prandtl 
nnnil'>er,  component  oiientation  and  chamber  width.  'I  his  coi  relation  was  found  to  be 
accniale  to  within  IDb  of  the  original  curve  fit  data.  Ifcat  transfer  in  I  (N75  was  found 
lo  occur  mainb  by  convection  arising  from  buoyancy  forces,  regardless  of  chamber 
uiilth.  Heat  transfer  in  r(N43  ua;  found  to  occur  mainly  by  molecular  diffusion  for 
chambci  wiilths  of  Ilium  or  greater,  and  by  convection  at  a  chamber  width  of  7mni. 
1  he  maximum  nnceitaipty  in  the  f3iissclr  and  Rayleigh  numbers  was  2..3'/'o,  based  on  a 
/froth  Older  unceifainty  antilysis,  and  occuied  at  the  lowest  power  level,  where  the 
maxiimnii  unceitainty  in  the  tempciatnie  measuiements  resided. 
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NOMENCLATURE 


Symbol  Description 
A  Area 

Total  wetted  surface  area 

Specific  heat 
5  Uncertainty 

emf  Thermocouple  voltage 

g  Gravitational  acceleration 

Gr  Grashof  number 

h  Average  heat  transfer  coefficient 

/y  Liquid  thermal  conductivity 

kp  Plexiglass  thermal  conductivity 

L  Component  length  in  the  vertical  direction 

NU[  Nusselt  number 

Pr  Prandtl  number 

Power  Power  supplied  to  a  heating  element 

Heat  loss  through  the  plexiglass  substrate 
Q„„  Net  power  dissipated  by  a  component 

Rc  Total  thermal  resistance  of  the  substrate 
Electrical  resistance  of  a  precision  resistor 
Ra^  Flux  based  Rayleigh  number 

ROf  Temperature  based  Rayleigh  number 

T  Temperature 

Average  component  temperature 
Tj  back  surface  temperature 

Tjy„  Film  temperature 

T„„i,  Average  heat  exchanger  temperature 
Volt  Input  voltage  from  the  DC  power  supply 
X  Non-dimensional  chamber  width 


Units 

nP 

nP 

J 

kg°C 

various 

volt 

m 

dimensionless 

IV 

m-°C 

IV 

nPC 

IV 

m°C 

m 

dimensionless 

dimensionless 

Watts 

Watts 

Watts 

°C 

IV 

Q 

dimensionless 

dimensionless 

°C 

“C 

’C 

“C 

“C 

volts 

dimensionless 


X 


Chamber  width 


Thermal  difTusivity 

Thermal  expansion  coelTicient 

Density 

Kinematic  viscosity 


I.  INTRODUCTION 


A.  PROBLEM  STATEMENT 

As  todays  circuit  technology  moves  beyond  Very  Large  Scale  Integration  (VLSI), 
increases  in  the  number  of  devices  per  chip,  package  densities  per  module,  and  power 
densities  at  both  the  chip  and  module  level  have  made  it  difficult  to  maintain  adequate 
operating  temperatures  (Incopera,  1990).  As  junction  temperature  increases,  microelec¬ 
tronic  chip  failure  increases  exponentially  (Park  and  Bergles,  1987).  Conversely,  for  ev¬ 
ery  20°C  drop  in  junction  temperarure,  chip  failure  rates  are  reduced  by  50% 
(Oktay,  1986).  Forced  air  cooling  and  indirect  liquid  cooling  appear  to  be  inadequate. 
More  emphasis  is  therefore  being  placed  on  direct  immersion  cooling  using  dielectric 
liquids.  This  may  involve  natural,  forced  or  mixed  convection  processes,  and  can  occur 
in  a  single  phase  as  well  as  pool  or  forced  boiling  convection  modes  (Incopera,  1990). 

Natural  convection  direct  immersion  cooling  using  dielectric  liquids  offers  a  high 
power  dissipation  capability,  while  providing  high  reliability  with  little  or  no  noise  (Liu, 
1987).  With  the  many  possible  applications  of  packaging  electronic  components  within 
an  enclosure,  few  studies  have  been  conducted  using  discrete,  protruding  heat  sources 
immersed  in  a  dielectric  liquid  (Incopera,  1988).  Studying  an  array  of  protruding  heat 
sources  of  uniform  size  and  simple  geometric  orientation  may  provide  valuable  insight 
and  prediction  capabilities  for  the  complex  geometric  arrangements  and  size  variations 
often  encountered  in  circuit  technology  today. 

B.  LITERATURE  REVIEW 

Park  and  Bergles  (1987)  conducted  natural  convection  heat  transfer  e.xperiments  in¬ 
volving  both  flush  and  protruding,  thin  foil  heaters  mounted  on  a  vertical  wall  within 
an  insulated  enclosure.  Using  water  and  R-113  as  working  fluids,  the  e.xperiment  con¬ 
sisted  of  three  different  cases: 

•  A  single  flush  mounted  heater 

•  An  array  of  flush  mounted  heaters  (in-line  and  staggered) 

•  An  array  of  protruding  heaters 

For  the  single  flush  mounted  heaters,  two  different  heights  (5mm  and  lOmm)  with 
variations  in  width  from  2mm  to  70mm  were  investigated.  The  height  and  width  of  the 
flush  mounted  heaters  were  limited  to  5mm  when  they  were  arranged  in  an  array.  The 


1 


protruding  heaters  were  4.9ntiii  in  the  vertical  direction,  5.3rnm  wide  and  extended 
1.1  nmi  from  the  substrate  surface.  The  experimental  results  showed  that  for  a  single 
flush  mounted  heater; 

•  The  heat  transfer  coefficient  increases  as  the  heater  width  decreases.  A  5mm  x 
5mm  test  section  displayed  a  heat  transfer  coefficient  80  to  100  times  larger  than 
the  widest  section  (70mm).  A  similar  eflect  was  observed  for  a  9.86mm  x  5mm  test 
section. 

•  The  width  effect  was  much  less  pronounced  for  w'ater  than  for  R-1 13. 

•  The  heat  transfer  coefficients  obtained  for  the  widest  test  section  were  20%  larger 
than  those  predicted  by  the  Fujii  and  Fujii  correlation  (Fujii  and  Fujii,  1976),  which 
is  based  on  boundary  layer  solutions  for  a  vertical  plate  with  constant  heat  flux. 

For  an  array  of  in-line  or  staggered  flush  mounted  heaters,  where  the  distance  betw'een 
them  was  varied,  the  results  showed  that: 

•  The  top  heater  had  a  substantially  lower  heat  transfer  coefficient  than  the  bottom 
heater. 

•  The  heat  transfer  coeficient  for  the  top  heater  increased  as  the  distance  between  the 
heaters  increased  up  to  a  value  of  G/H  =  3.5,  wirere  H  is  the  height  of  the  heater 
and  G  the  distance  between  them.  For  values  of  G/H  greater  than  3.5,  the  heat 
transfer  coefficient  ratio  remained  constant. 

•  For  staggered  flush  mounted  heaters,  the  heat  transfer  coefficient  increased  as  the 
transverse  distance  between  the  heaters  increased  at  small  G/TI.  The  opposite  ef¬ 
fect  occured  at  large  G/H. 

For  a  single  protruding  heater,  the  heat  transfer  coefficient  was  found  to  be  14°/o 
larger  than  the  flush  mounted  value,  and  120%  larger  than  the  predicted  values  using  the 
Fujii  and  Fujii  correlation.  Heat  transfer  coefficients  for  the  top  heater  in  an  array  were 
higher  than  the  bottom  and  increased  as  the  distance  between  them  increased.  This  w'as 
in  contrast  to  the  result  obtained  for  flush  mounted  heaters. 

Kelleher  et  al.  (1987)  conducted  an  e.xperimental  investigation  of  natural  convection 
in  a  W'ater  filled,  rectangular  enclosure  with  a  single  heated  protrusion  on  a  vertical  wall. 
All  vertical  exterior  surfaces  W'ere  insulated  and  assumed  adiabatic.  The  upper  and  lower 
horizontal  boundaries  w'ere  maintained  at  an  isothermal  sink  temperature,  .  The 
heated  protrusion  had  the  following  dimensions:  H  =  25.4mm  W=  203.0mm  and 
protruded  12.8mm  from  the  surface  of  the  wall.  The  enclosure  w'as  constructed  such  that 
the  heated  protrusion  could  be  set  at  3  different  elevations  within  the  enclosure.  Ex¬ 
perimental  runs  W'ere  conducted  in  the  Rayleigh  number  range: 
.99  X  10‘  <  7?a  <  8  X  10^  A  flow  visualization  study  was  also  carried  out.  The  exper¬ 
imental  results  showed  that; 
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•  For  a  given  Rayleigh  number,  Nusseli  number  decreased  as  the  height  of  the  heater 
within  the  enclosure  increased. 

•  Flow  visualization  revealed  a  dual  celled  flow  which  consisted  of  a  buoyancy  driven 
upper  cell  and  a  shear  driven  lower  cell.  The  motion  of  the  lower  cell  arises  due  to 
the  viscous  drag  of  the  upper  cell.  The  majority  of  fluid  motion  takes  place  in  the 
upper  cell  and  accounts  for  the  most  of  the  convective  heat  transfer. 

Lee  et  al.  (1987)  conducted  a  numerical  simulation  of  laminar  natural  convection  in 
a  water-filled  rectangular  enclosure,  using  a  single  heated  protrusion  on  a  vertical  wall. 
This  was  a  companion  paper  to  the  research  done  by  Kelleher  et  al.  (1987).  A 
2-dimensional,  full  variable-property,  finite  difference  method  calculation  was  used  to 
simulate  the  experimental  conditions  investigated  by  Kelleher  et  al.  (1987).  The  numer¬ 
ical  solution  to  the  governing  conservation  equations  revealed  the  following: 

•  At  high  Rayleigh  numbers,  there  is  reasonable  agreement  between  the  experimental 
and  numerical  data  for  the  Nusselt  number. 

•  At  low  Rayleigh  numbers,  the  numerical  solution  predicts  values  of  the  Nusselt 
number  as  much  as  50%  higher  than  those  observed  experimentally. 

•  In  contrast  to  the  experimetal  results  of  Kelleher  et  al.  (1987),  the  numerical  sol¬ 
ution  showed  that  the  Nusselt  number  increased  as  the  elevation  of  the  heater 
within  the  enclosure  was  increased. 

•  Flow  visualization  showed  that  the  flow  is  concentrated  in  the  region  above  the 
heater.  The  region  below  the  heater  was  essentially  stratified.  Flow  in  the  upper 
region  was  buoyancy  driven,  while  the  region  below  the  heater  v'as  shear  driven 
by  the  buoyant  flow  above. 

Chen  and  Kuo  (1988)  studied  the  natural  convective  heat  transfer  from  10mm  x 
20mm  x  150mm  block-like  heat  generating  modules.  A  numerical  study  was  first  con¬ 
ducted  using  a  laminar,  2-dimensional,  constant  property  finite -difference  scheme. 
Reslults  were  provided  for  a  Raleigh  number  range  of:  10^  <  Ra<  10^  at  Pr  =  0.7  .  The 
numerical  results  were  corroborated  experimentally  by  measuring  temperatures  with  a 
laser  holographic  interferometer.  The  results  indicate  that: 

•  Heat  transfer  on  the  upstream  side  of  the  heated  surface  is  in  general  higher  than 
the  downstream  side. 

•  The  heat  transfer  coefficient  for  the  upper  heater  is  affected  by  the  flow  generated 
by  the  lower  heater.  This  effect  increases  as  the  Rayleigh  number  increases. 

Keyhani  et  al.  (1991)  conducted  experiments  to  determine  the  effect  of  aspect  ratio 
on  the  natural  convection  in  an  enclosure  with  5  protruding  heat  sources.  The  enclosure 
had  inner  dimensions:  H=  171.35mm  L=  141.0mm  W=  139.7mm,  and  was  constructed 
of  25.4mm  thick  plexiglass.  The  upper  horizontal  boundary  of  the  enclosure  was  main- 
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tained  at  a  contant  temperature.  The  width  of  the  enclosure  was  varied  from  13.5mm 
55mm  by  inserting  a  movable,  vertical  plate  into  the  enclosure.  The  protruding  heaters 
were  15mm  in  the  vertical  direction,  and  extended  9mm  from  the  surface  of  the  phenolite 
wall  on  which  they  were  mounted.  The  protruding  heaters  were  mounted  such  that 
identical  rows  of  protruding  heaters  and  flush  unheated  sections  were  obtained.  Using 
ethylene  glycol  as  a  working  fluid,  experimental  runs  were  carried  out  for  power  inputs 
of  2W  to  13  W  per  heated  section,  for  6  enclosure  widths.  The  Prandtl  number  for 
ethylene  glycol  varied  from  62  to  110.  Flow  visualization  experiments  were  also  con¬ 
ducted  to  investigate  the  flow  structure.  A  correlation  which  relates  the  local  Nusselt 
number  to  the  local  modified  Rayleigh  number  and  aspect  ratio  (~j^)  ''’as  reported: 

Nu^  =  0.296(Rayf^^\ 

where  Ra^  =  Ra^'Su.  This  correlation  is  independent  of  the  number  of  heaters  in  the 
array,  enclosure  width  to  protrusion  heater  height  ratio,  and  vertical  height  location  of 
the  heaters.  The  conditions  for  which  this  conclusion  is  valid  are  outlined  in  Keyhani 
et  al.  (1991).  The  averaged  deviation  of  this  correlation  from  the  experimental  data  was 
reported  as  4.5%. 

The  experimental  study  reported  in  this  thesis  is  a  direct  continuation  of  the  work 
conducted  by  Aytar  (1991),  who  investigated  the  heat  transfer  characteristics  of  a  3x3 
array  of  horizontally  oriented  components.  Frequent  references  are  made  to  Aytar's 
(1991)  research  throughout  this  study,  since  most  facets  of  his  study  have  been  repeated 
herein. 

Aluminum  blocks  24mm  x  8mm  x  6nim,  each  simulating  a  20  pin  dual-inline  pack¬ 
age,  were  mounted  to  a  plexiglass  substrate  to  produce  a  3x3  array  of  vertically  oriented 
simulated  electronic  components  (Figure  5  and  Figure  6).  The  resulting  circuit  board 
was  afiixed  in  a  rectangular  plexiglass  enclosure  and  immersed  in  a  dielectric  liquid.  The 
enclosure  had  inner  dimensions  of:  H=  152.0mm  L=  203.2mm  W=  82.6mm,  and  was 
constructed  of  25.4mm  thick  plexiglass  (Figure  1).  A  constant  heat  flux  was  supplied 
to  each  of  the  9  components  using  an  Inconel  foil  resistive  heating  element  connected 
in  parallel  to  a  0-100  0-5A  DC  power  supply.  Power  supplied  to  the  heating  elements 
varied  from  0.115VV  to  2.90VV  per  element.  The  effect  of  chamber  width  w’as  investigated 
at  the  following  values:  42,  30,  18,  11  and  7mm.  The  upper  horizontal  boundary  of  the 
enclosure  was  maintained  at  lO^C  using  a  circulating  chilled  water  heat  exchanger.  All 
other  exterior  surfaces  were  insulated.  The  distance  from  the  front  center  face  of  the 


4 


difTerent  thicknesses.  Temperature  measurements  at  the  component  faces,  the  rear  plane 
of  the  the  circuit  board,  and  the  heat  exchanger  boundar\’,  were  made  using  copper- 
constantan  thermocouples.  The  thermocouples  were  tack  welded  in  small  wells  below 
the  aluminum  block  surface  to  minimize  any  contact  resistance  (Figure  3  and 
Figure  4).  An  HP3497A  data  acquisition  system,  controlled  by  an  HP9826  microcom¬ 
puter,  was  used  to  collect  the  raw  temperature  data  after  the  system  achieved  steady 
state.  Steady  state  conditions  were  assumed  when  surface  temperature  fluctuations  were 
less  than  1°C  in  10  minutes  for  FC-75  and  0.4°C  for  FC-43.  These  criterion  were  based 
on  steady  state  surface  temperature  fluctuation  measurements  made  by  Aytar  (1991)  for 
a  horizontal  component  orientation. 

C,  OBJECTIVES 

The  objectives  of  this  work  were: 

1.  Manufacture  a  3x3  array  of  vertically  oriented  simulated  electronic  components. 

2.  Collect  data  for  the  vertical  component  orientation  by  varving  the  power  level 
supplied  to  the  heating  elements  and  the  width  of  the  enclosure  using  FC-75  and 
FC-43  as  working  fluids. 

3.  Reduce  the  data  into  both  dimensional  and  non-dimensional  heat  transfer  param¬ 
eters  for  the  chamber  widths  and  power  levels  investigated. 

4.  Develope  an  emperical  correlation  which  predicts  the  heat  transfer  characteristics 
of  the  array,  taking  into  account  the  effects  of  Rayleigh  number,  Prandtl  number, 
chamber  width  and  component  orientation. 

5.  Make  comparisons  between  the  vertical  component  orientation  and  the  horizontal 
component  orientation,  using  both  dimensional  and  non-dimensional  data. 

6.  Based  on  the  findings  of  this  study,  make  recommendations  for  future  research. 
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II.  EXPERIMENTAL  APPARATUS 


A.  TEST  CHAMBER 

During  the  course  of  this  experimental  study,  2  test  chambers  were  used.  Each  was 
constructed  of  plexiglass  with  inner  dimensions;  L=  203.2mm  W=  82.6mm 
H=  152.0mm,  and  a  12mm  thick  bottom  boundary  (Figure  1).  The  first  chamber  had 
a  wall  thickness  of  12.7mm,  and  was  used  to  obtain  data  for  FC-75  at  chamber  widths 
of  42,  30,  18  and  1 1mm.  During  these  experimental  runs  the  chamber  experienced  brittle 
fractures  at  numerous  points  along  its  boundary,  particularly  at  the  corners.  This  re¬ 
sulted  in  the  chamber  having  to  be  discarded.  A  photograph  of  the  first  enclosure,  with 
the  chilled  water  heat  exchanger  mounted  in  place,  is  presented  in  Figure  2  to  provide 
some  general  insight  into  the  enclosure's  detail.  A  second  chamber  was  constructed  with 
the  following  modifications; 

•  The  wall  thickness  was  increased  to  25.4mm 

•  The  lower  steel  tubing,  which  served  as  a  fill  and  removal  point  for  the  dielectric 
liquid,  was  moved  slightly  forward  of  the  circuit  board  surface,  to  provide  an  un¬ 
obstructed  flow  path  for  the  dielectric. 

•  The  upper  steel  tubing,  which  cannot  be  seen  if  Figure  2,  was  removed  since  it  no 
longer  served  any  useful  purpose.  It  was  previously  used  as  an  air  removal  point 
in  a  past  enclosure  design. 

•  Threaded  studs  were  permanently  set  within  the  chamber  walls  in  order  to  lessen 
the  wear  on  the  threads  within  the  plexiglass  walls  themselves.  The  heat  exchanger 
was  then  fastened  to  the  chamber  using  washers  and  nuts. 

A  top  view  of  the  second  chamber  with  the  circuit  board  and  vertically  mounted  com¬ 
ponents  are  shown  in  Figure  1.  The  experimental  runs  carried  out  for  FC-75  using  the 
initial  chamber  were  also  completed  using  the  second  chamber. 

B.  HEAT  EXCHANGER 

The  upper  horizontal  boundary  of  the  chamber  was  maintained  at  lO^C  using  a 
water-cooled  channel  type  heat  exchanger.  The  shell  was  constructed  of  plexiglass  and 
contained  a  3mm  thick  aluminum  plate,  which  provided  a  nearly  uniform  temperature 
surface  between  the  dielectric  liquid  and  the  cooling  water.  A  detailed  drawing  and  de¬ 
scription  of  the  heat  exchanger  and  its  construction  was  given  by  Torres  (1988).  Prior 
to  this  study,  the  aluminum  plate  contained  2  rows  of  3  symmetrically  placed 
thermocouples.  An  additional  row  of  thermocouples,  as  recommended  by  Aytar  (1991), 
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was  added  to  the  aluminum  plate  to  provide  a  better  measurement  of  the  average  upper 
boundry  temperature  at  the  various  chamber  widths.  Flow  through  the  channels,  as  well 
as  the  inlet  water  temperature,  was  adjusted  from  the  refrigerated  circulating  bath  unit. 

C.  SIMULATED  ELECTRONIC  COMPONENTS 

Rectangular  aluminum  blocks  24mm  x  8mm  x  6mm  were  used  to  simulate  20  pin 
dual-in-line  packages.  Prefabricated  copper-constantan  thermocouples  were  used  to 
measure  the  surface  temperature  of  the  components.  A  detailed  drawing  of  the  alumi¬ 
num  blocks  and  the  thermocouple  wells  is  shown  in  Figure  3  (after  Benedict,  1988).  The 
thermocouples  were  mounted  to  the  aluminum  blocks  in  accordance  with  the  procedure 
described  by  Aytar  (1991).  A  photograph  of  a  finished  component  is  provided  in 
Figure  4. 

D.  HEATING  ELEMENTS 

Inconel  foil  resistive  heating  elements  were  used  to  provide  a  constant  heat  flux  to 
the  back  surface  of  the  components,  via  a  O-IOOV  0-5A  DC  power  supply.  The  heaters 
were  mounted  to  the  plexiglass  substrate  in  accordance  with  the  procedure  described  by 
Aytar  (1991).  The  resistance  of  each  heating  element  was  approximately  IIH.  Each 
heating  element  was  connected  in  series  with  a.  2Q.  ±  2.5%  precision  resistor.  The  9 
heater-precision  resistor  series  elements  were  then  connected  in  a  parallel  combination 
with  the  DC  power  supply.  By  measuring  the  voltage  drop  across  each  component  for 
a  given  power  supply  voltage,  the  current,  and  therefore  the  power  provided  to  each 
component,  was  determined. 

E.  CIRCUIT  BOARD 

The  aluminum  blocks  were  mounted  to  a  plexiglass  substrate  with  their  largest  ver¬ 
tical  dimension  (24mm)  oriented  in  the  direction  of  the  gravitational  vector,  to  produce 
a  3x3  array  of  simulated  electronic  components.  The  substrate  was  152.0mm  in  the 
vertical  direction,  203.2mm  wide  and  12mm  thick.  Side  walls  were  added  to  the  back  of 
the  circuit  board  to  provide  stability  when  it  was  inserted  into  the  test  chamber.  The 
dimensions  of  the  side  walls  were  such  that  when  the  board  was  inserted  into  the  cham¬ 
ber  without  a  spacer,  the  distance  from  the  surface  of  the  board  containing  the  compo¬ 
nents,  to  the  inner  surface  of  the  chamber,  was  42mm.  A  small  groove  was  made  in  the 
upper  right  hand  comer  of  the  board  to  allow  air  to  pass  from  the  front  of  the  chamber 
to  the  back,  where  it  could  then  be  vented  to  the  atmosphere  as  the  chamber  was  filled. 
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A  front  and  plan  view  of  the  circuit  board  with  the  vertically  mounted  components  are 
shown  in  Figure  5.  A  photo  of  the  finished  circuit  board  is  shown  in  Figure  6. 

F.  SYSTEM  HARDWARE 


All  thermocouples  and  heating  elements  were  connected  to  an  HP3497A  data  ac¬ 
quisition  system,  which  was  controlled  by  an  HP9826  microcomputer.  The  channel 
numbering  scheme  was  the  same  as  used  by  Aytar  (1991),  and  is  repeated  below: 

•  Channels  0  through  53:  Component  surface  temperatures 

•  Channels  54  through  56:  Heat  exchanger  surface  temperatures 

•  Channels  57  through  60:  Circuit  board  back  plane  temperatures 

•  Channels  61  through  70:  DC  power  supply  and  all  heating  elements 

•  Channels  70  through  75:  Circuit  board  back  plane  temperatures 

•  Channel  76:  Ambient  temperature 


Components  in  the  vertical  orientation  were  numbered  in  the  same  way  as  in  the  hori¬ 
zontal  orientation:  bottom  to  top  and  right  to  left,  viewed  from  the  surface  containing 
the  components.  The  horizontal  orientation  was  rotated  90°  to  arrive  at  the  vertical 
orientation.  The  component  faces  were  then  numbered  using  the  scheme  described  by 
Aytar  (1991):  center,  top,  right,  left,  bottom  and  back,  beginning  with  component  #l  and 
ending  wth  component  #9. 
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Figure  1.  Top  View  of  the  Test  Chamber  and  Circuit  Board  with  Vertically  Oriented 
Components. 
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Figure  2.  Plexiglass  Enclosure  with  the  Water  Chilled  Heat  Exchanger  Mounted 
in  Place. 
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All  dimensions  in  millimeters! 


Figure  3.  Aluminum  Blocks  and  Thermocouple  Wells. 
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Figure  4.  Aluminum  Block  Mounted  with  Copper-Constanlan  Thermocouples  and 
an  Inconel  Foil  Heating  Element. 
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Figure  5.  Front  and  Plan  View  of  Circuit  Board  and  Vertically  Oriented  Compo¬ 
nents. 
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Figure  6.  Circuit  Board  with  Vertically  Oriented  Components. 
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III.  EXPERIMENTAL  PROCEDURE  AND  DATA  ANALYSIS 


A.  SYSTEM  PREPARATION 

Prior  to  each  experimental  run,  the  following  preparatory  steps  were  completed. 

1.  The  required  spacer  (if  any)  was  inserted  into  the  chamber.  The  heat  exchanger 
was  affixed  to  the  chamber  only  enough  to  ensure  that  the  circuit  board  and  spacer 
were  secure.  Using  a  ruler  graduated  in  millimeters,  the  chamber  width  was  meas¬ 
ured  at  various  points  along  the  height  of  the  plexiglass  substrate  from  the  exterior 
of  the  chamber,  to  ensure  that  the  distance  from  the  circuit  board  to  the  spacer 
was  consistent.  These  measurements  were  taken  on  both  sides  of  the  chamber. 
Once  the  chamber  width  was  set  at  the  proper  distance,  the  heat  exchanger  was 
then  completely  affixed  to  the  chamber. 

2.  The  heat  exchanger  supply  and  return  lines  were  attached,  and  the  system  tested 
for  any  interior  or  exterior  leakage.  This  test  was  conducted  with  the  How  setting 
at  its  maximum. 

3.  Upon  a  satisfactory  heat  exchanger  seal  test,  the  flow  control  was  set  to  zero  and 
the  bath  temperature  control  set  to  lO'C.  The  bath  was  allowed  to  reach  test 
conditions  while  the  remaining  preparations  were  completed. 

4.  The  chamber  was  filled  with  the  appropriate  dielectric  liquid,  ensuring  that  as  much 
air  as  possible  was  removed  from  the  chamber  (this  is  especially  important  at 
smaller  chamber  widths).  The  chamber  was  observed  for  any  visible  sign  of 
leakage. 

5.  Foam  insulation  padding  was  affixed  to  the  bottom  and  vertical  surfaces  of  the 
chamber. 

6.  The  appropriate  heat  exchanger  thermocouples  were  chosen  based  on  the  chamber 
width  under  investigation. 

7.  All  thermocouples  and  heating  elements  were  checked  for  continuity  by  stepping 
through  each  channel  of  the  HP3497A  data  acquisition  system  and  monitoring  the 
digital  voltmeter. 

B.  EXPERIMENTAL  PROCEDURE 

Once  system  preparations  were  complete,  the  following  steps  were  followed  to  ac¬ 
complish  an  experimental  run,  which  consisted  of  a  single  chamber  width  at  all  power 
levels. 

1.  The  voltage  on  the  DC  power  supply  was  set  to  produce  the  required  power  level. 

2.  The  refrigerated  bath's  flow  control  was  opened  to  its  maximum  position. 

3.  The  heat  exchanger  thermocouples  (channels  54-56)  were  monitored,  and  the  bath 
temperature  control  adjusted  to  attain  the  10®C  upper  boundary’  condition. 

4.  Once  the  average  heat  exchanger  temperature  stabilized,  an  initial  data  set  was 
taken  to  establish  a  system  start  time. 
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5.  Steady  state  conditions  were  assumed  when  the  temperature  fluctuations  from  all 
thermocouples  were  less  than  TC  in  10  minutes  for  FC-75  and  0.4'’C  in  10  minutes 
for  FC-43.  These  conditions  were  the  same  as  those  used  by  Aytar  (1991),  and  are 
based  on  steady  state  surface  temperature  fluctuation  measurements.  As  reported 
by  Aytar  (1991),  the  time  constant  of  the  system  was  a  function  of  chamber  width 
and  the  dielectric  fluid  used.  This  also  held  true  for  the  vertical  component  orien¬ 
tation.  Typical  time  requirements  to  achieve  steady  state  from  ambient  conditions 
were  approximately  3  1/2-4  hours.  For  subsequent  runs  at  the  same  chamber 
width,  the  time  decreased  to  approximately  2-2  1/2  hours. 

C.  DATA  ACQUISITION 

Data  acquisition  and  processing  was  accomplished  using  the  software  programs 
ACQUIRE  and  CALCDIEL.  These  programs  were  the  same  as  those  used  by  Pamuk 
(1987),  Benedict  (1988),  Torres  (1988),  Powell  (1989)  and  Aytar  (1991),  with  the  follow¬ 
ing  modifications: 

•  A  one  point  calibration  was  done  for  each  card  of  the  data  acquisition  system  using 
an  ice  bath  reference.  Three  random  channels  per  card  were  tested  against  the  ice 
bath  reference,  and  the  arithmatic  average  of  the  3  cmf  readings  was  assumed  to 
be  a  constant  offset  for  that  particular  data  acquisition  system  card. 

•  CALCDIEL  was  made  completely  general  and  interactive,  so  that  any  combination 
of  power  level,  chamber  width,  dielectric  liquid  or  component  orientation  could  be 
used  as  an  input. 
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IV.  DATA  ANALYSIS 


The  following  analysis  was  used  to  reduce  the  raw  data  obtained  by  ACQUIRE  into 
non-dimensional  parameters  (Nu,,  Ra„  Ra^.  These  non-dimensional  parameters  were 
determined  for  each  component  at  the  various  power  levels  and  chamber  widths  for  both 
FC-75  and  FC-43  using  CALCDIEL.  The  assumptions  used  to  arrive  at  the  expressions 
in  the  following  analysis  are  given  below: 

•  Heat  transfer  from  the  heating  elements  to  the  plexiglass  substrate  occured  via 
1 -Dimensional  conduction. 

•  A  linear  temperature  profile  existed  across  the  thickness  of  the  plexiglas  substrate. 

•  The  surface  area  of  the  plexiglass  substrate  not  occupied  by  a  component  was 
adiabatic  (i.e.  there  was  no  heat  transfer  between  the  dielectric  liquid  and  the 
substrate). 

•  The  thermophysical  properties  of  the  dielectric  liquid  near  each  component  were 
uniform  and  constant,  and  were  evaluated  at  the  film  temperature 

•  The  thermal  conductivity  of  the  plexiglass  substrate  was  homogenous. 

•  There  was  negligable  contact  resistance  between  the  aluminum  blocks  and  the  foil 
heaters. 

•  There  was  negligable  contact  resistance  between  the  foil  heaters  and  the  plexiglass 
substrate. 

A  components  av'erage  surface  temperature  was  approximated  as  an  area  based  av¬ 
erage  over  the  entire  wetted  surface  area  of  the  component: 


s 


where, 

A„T,  =  area  and  temperature  of  the  /'*  surface 
==  total  wetted  surface  area  of  a  component 

The  net  power  dissipated  by  a  heating  element  was  determined  by  measuring  the 
voltage  drop  across  the  heating  element,  and  multiplying  it  by  the  current  passing 
through  the  precision  resistor-heating  element  series  combination: 
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where, 

emf  =  voltage  drop  across  a  resistive  heating  element 

volt=  DC  power  supply  voltage 

Rp  =  resistance  of  a  precision  resistor  (2Q) 

The  total  heat  loss  from  the  array's  heating  elements  to  the  plexiglass  substrate  was 
calculated  as  a  sum  of  contributions  from  each  independent  heating  element.  The 
process  was  considered  to  occur  via  1-D  conduction  across  the  thickness  of  the 
substrate,  where  a  linear  temperature  profile  was  assumed. 
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Qloss 


AT=T,,-T,p 


where, 

=  component  back  surface  temperature 
=  the  back  circuit  board  temperature 
=  conductive  thermal  resistance  of  the  plexiglass 
=  thermal  conductivity  of  the  plexiglass 
A  =  rear  face  surface  area  of  the  comonent 
t  =  thickness  of  the  circuit  board 

The  net  heat  transfered  from  a  component  to  the  surrounding  dielectric  liquid  was 
determined  by  subtracting  the  total  heat  loss  to  the  substrate  via  1-D  conduction  from 
the  net  power  dissipated  by  its  heating  element: 


=  POWER 
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The  average  heat  transfer  coefficient  for  a  component,  h,  was  determined  from 
Newton's  Law  of  Cooling  for  convective  heat  transfer; 


Q. 


net 


^loia^'t'avg  ^sink) 


Tsink  — 


^54  +  ^’55  +  r. 


56 


3 


where, 

T'tink  ~  average  temperature  of  the  heat  exchanger 
The  non-dimensional  Nusselt  number  (Nu,),  was  defined  as; 


V  hL 

Nu,=  — 


where, 

L  =  vertical  dimension  of  the  component  in  the  direction  of  the  gravitational  vector 
kf  =  thermal  conductivity  of  the  dielectric  liquid 

The  non-dimensional  Grashof  number,  which  relates  the  effects  of  buoyant  forces 
to  viscous  forces  in  a  natural  convection  environment,  was  defined  as; 

G,  = - ^ - 

The  thermal  diffusivity,  a,  which  indicates  the  thermal  diffusion  characteristics  of  a 
fluid,  was  found  from; 


a  = 


pCp 


The  Prandtl  Number,  which  relates  the  relative  rates  of  viscous  to  thermal  diffusion 
for  a  fluid,  was  defined  as: 
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The  temperature  based  Rayleigh  number,  which  is  the  product  of  the  non- 
dimensional  Grashof  and  Prandtl  numbers,  is  synonymous  with  the  Reynolds  number  ( 
Re)  used  in  forced  convection  heat  transfer.  This  parameter  was  defined  as; 


ROf  =  Gr »  Pr  = 


va 


AT—  Tayg  —  Tjtnk 

The  flux  based  Rayleigh  number  was  defined  as: 


The  thermophysical  properties  of  the  dielectric  liquids  were  determined  from  the 
following  relations,  and  evaluated  at  the  film  temperature, 


'Tfilm  — 


Tavg  T  T^slnk 


Thermal  Conductivity,  k. 


for  FC-75: 


for  FC-43:  kf  = 


[0.65  -  7.8947  x  10"^  x  Tyj/J 


10 


[0.666  -  9.864  x  10"^  x 


10 


Density,  p 


for  FC-75:  p  =  [1.825  -  0.00246  x  7}5/J  x  10^ 
for  FC-43;  p  =  [1.913  -  0.00218  x  7}j/J  x  10^ 


Specific  Heat,  C 


for  FC-75  and  FC-43:  =  [0.241 1 1 1  +  3.7037  x  10"^  x  7};/J  x  4.18  x  10^ 
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Kinematic  Viscosity,  v 
FC.75: 

v  =  [l. 4074  -2.964  x  10"^  x  + 3.8018  x  10“*  x  7^, 2.7308  x  l0“‘x  7^,„  + 8.1679  x  10“’  x  7^,„]  x  10“ 


FC-43: 

v=  [8.8750  -0,47007  x  ?},/„+  1.387  x  10“’ x  7^;„- 2.1469  x  10“*  x  +  1.3139  x  10“‘ x  x  10“‘ 


Thermal  Expansion  Coefficient,  p 

for  FC-75:  /?  =  [ 

for  FC-43:  /?  =  [ 


0.00246 


(1.825  - 0.00246  X  Tyj!„) 

_ 0.00218 _ 

(1.913  -0.00218  X  TyjyJ 
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V.  RESULTS 


A.  GENERAL 

Dimensional  and  non-dimensional  data  output  was  used  to  study  the  heat  transfer 
characteristics  of  a  3x3  array  of  vertically  oriented  simulated  electronic  components  us¬ 
ing  two  dielectric  fluids:  FC-75  and  FC-43,  at  chamber  widths  of  42,  30,  18,  11,  and  7 
mm  and  power  levels  from  0.1 14-2. 9W.  The  non-dimensional  data  from  this  study  was 
combined  with  data  obtained  by  Aytar  (1991),  to  produce  an  empirical  correlation  which 
accounted  for  variations  in: 

•  Rayleigh  number 

•  Chamber  Width 

•  Prandtl  number 

•  Component  orientation 

Each  of  these  effects  were  addressed  separately  and  then  subsequently  combined  to 
develope  a  relation  of  the  form: 

Nu,  = 

where  Nui  is  the  Xusselt  number  based  on  the  components  dimension  in  the  direction 
of  the  gravitational  vector,  X  is  a  non-dimesional  chamber  width  based  on  the  maximum 
chamber  width  attainable  and  a,  b^,  b-^  and  b-^  are  constants. 

B.  DIMENSIONAL  RESULTS 

Figure  7  and  Figure  8  show  the  average  array  temperature  as  a  function  of  both  the 
net  power  dissipated  by  the  array  and  chamber  width,  for  both  the  vertical  and  hori¬ 
zontal  orientations.  For  both  orientations,  average  array  temperatures  are  seen  to  be 
independent  of  chamber  width  at  power  levels  less  than  O.IW,  and  nearly  independent 
of  chamber  width  at  all  power  levels  for  chamber  widths  of  30mm  or  greater.  At  power 
levels  greater  than  0.5W,  the  data  form  straight  lines  with  slopes  that  increase  as  the  wall 
spacing  decreases  from  42mm.  This  trend  was  also  observed  for  FC-43.  As  mentioned 
previously,  two  data  sets  were  taken  for  FC-75  at  chamber  widths  of  42,  30,  18  and 
1 1mm.  When  the  data  for  the  2  test  chambers  were  plotted  as  in  Figure  7,  the  two  data 
sets  fell  within  the  maximum  uncertainty  range  of  the  thermocouples.  The  uncertainty 
was  specified  as  1%  of  the  temperature  value  or  ±0.5‘’C,  whichever  was  greater.  The 
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increase  in  wall  thickness  therefore  had  no  effect  on  the  data.  The  two  data  sets  were 
averaged  to  produce  one  data  set  for  those  chamber  widths.  The  averaged  data  set  was 
used  to  produce  all  quantitative  results,  whereas  all  qualitative  results  were  based  on  the 
data  taken  using  the  second  chamber. 

The  ma.ximum  increase  in  the  average  array  temperature  for  the  vertical  orientation, 
for  all  chamber  wdths  and  power  levels,  was  34.9'’C  for  FC-75  and  dS.S'C  for  FC-43. 
For  the  horizontal  orientation,  the  maximum  temperature  increase  was  42.5°C  for  FC-75 
and  44.5®C  for  FC-43.  Maximum  average  component  temperatures  for  the  vertical  ori¬ 
entation  were  52°C  for  FC-75  and  68°C  for  FC-43,  whereas  for  the  horizontal  orien¬ 
tation  maximum  average  component  temperatures  were  59.8°C  for  FC-75  and  68.3‘’C  for 
FC-43.  These  maximum  temperatures  occurred  at  a  power  level  of  approximately  2.8W 
and  a  chamber  width  of  7mm  for  both  orientations  and  dielectric  liquids. 

The  average  component  temperatures  for  FC-75  and  a  vertical  component  orien¬ 
tation  were  nearly  uniform  throughout  the  array  for  power  levels  less  than  1.70W  and 
chamber  widths  of  42,  30,  and  18mm.  As  the  chamber  width  was  decreased  to  11mm, 
the  power  level  below  which  component  average  temperatures  remained  uniform,  de¬ 
creased  to  0.80W.  At  7mm,  this  power  level  decreased  to  0.34W.  For  each  chamber 
width,  temperature  differences  between  rows  developed  above  those  power  levels  where 
the  average  component  temperatures  remained  uniform.  These  temperature  differences 
became  more  pronounced  as  the  power  level  increased,  with  noticable  patterns  forming 
between  the  rows.  However,  average  component  temperatures  within  each  row  re¬ 
mained  nearly  uniform.  In  order  to  make  comparisons  between  the  rows,  an  average 
row  temperature  was  determined  from  the  3  components  within  each  row.  For  FC-75 
and  a  vertical  component  orientation,  the  order  of  average  row  temperatures  was: 
middle  >  top  >  bottom.  This  pattern  occured  at  all  chamber  widths,  but  at  the  different 
power  levels  where  the  average  component  temperatures  ceased  to  be  uniform.  This 
pattern  may  be  attributed  to  the  dominance  of  the  buoyancy  forces  over  viscous  forces. 
A  strong  upward  flow  induced  by  the  buoyancy  forces  retards  the  developement  of  the 
ascending  boundary  layer.  The  decrease  in  the  average  row  temperature  from  the  middle 
to  top  row  was  attributed  to  the  close  proximity  of  the  components  in  the  top  row  to 
the  upper  isothermal  boundary.  The  maximum  temperature  difference  between  the 
middle  and  bottom  rows  was  6.8'’C,  and  occured  at  a  power  level  of  2.86W  and  a 
chamber  width  of  1 1mm.  There  was  no  component,  at  any  power  level  or  wall  spacing, 
which  was  consistently  at  either  the  maximum  or  minimum  average  temperature.  The 
only  conclusion  that  could  be  made  was  that  the  maximum  temperature  occured  on  a 
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component  in  the  middle  row,  and  the  minimum  on  a  component  in  either  the  top  or 
bottom  row. 

For  FC-43  and  a  vertical  component  orientation,  average  component  temperatures 
remained  nearly  uniform  for  all  chamber  widths  at  power  levels  of  O.S^W  or  less.  Again, 
temperature  differences  between  rows  developed  at  power  levels  greater  than  0.34W, 
with  average  component  temperatures  within  each  row  remaining  nearly  uniform.  The 
order  of  average  row  temperatures  for  FC-43  and  a  vertical  component  orientation  was: 
top  >  middle  >  bottom.  This  pattern  occured  at  all  chamber  widths  and  power  levels, 
with  the  exception  of  a  7mm  chamber  width  and  power  levels  of  2.24W  and  2.86W, 
w’herc  it  was:  middle  >  top  >  bottom.  For  the  pattern  mentioned  first,  viscous  forces 
dominate  buoyancy  forces  and  the  circuit  board  acts  like  a  vertical,  flat  plate  with  a 
constant,  uniform  heat  flux  (Blasius  flow).  Here,  the  ma.ximum  temperature  difference 
between  the  top  and  bottom  average  row  temperatures  w'as  3.5“C,  and  occured  at  a 
chamber  width  of  7mm  and  a  power  level  of  1.70W.  The  maximum  average  component 
temperature  in  the  array  occured  almost  exclusively  at  chip  6,  and  was  always  in  the 
upper  row.  For  the  second  pattern  mentioned  above,  which  w^as  the  same  as  observed 
for  FC-75  and  a  vertical  component  orientation,  buoyancy  forces  appear  to  have  over¬ 
come  viscous  forces.  The  relatively  small  distance  from  the  center  front  face  of  the 
components  to  the  wall  (1mm)  may  have  also  affected  the  developement  of  the  ascend¬ 
ing  boundary'  layer,  as  well  as  the  upward  flow  induced  by  the  buoyancy  forces.  The 
maximum  temperature  difference  between  the  middle  and  bottom  row  average  temper¬ 
atures  was  6.4'’C  and  occured  at  a  power  level  of  2.86W.  The  maximum  average  com¬ 
ponent  temperatures  for  these  tw'o  power  levels  at  a  7mm  chamber  width,  occured  at 
chip  2.  The  minimum  average  component  temperature  occured  at  chip  7  for  all  cases 
described  above. 

C.  NON-DIMENSIONAL  RESULTS 
1.  FC-75 

The  maximum  deviation  in  the  flux  based  Rayleigh  (Ra/)  number  and  Nusselt 
number  (Nu,),  for  a  vertical  component  orientation,  w'-as  9.6%  and  18.8%  respectively. 
Both  occured  at  a  chamber  width  of  7mm  and  a  power  setting  of  2.87\V.  The  average 
deviation  in  Ra^  and  Nu,,  among  all  the  data  points  for  the  vertical  orientation,  was 
4.5%  and  10.8%  respectively.  For  the  horizontal  component  orientation,  Aytar  (1991) 
reported  a  maximum  deviation  of  7%  in  both  Ruyand  Nu,  for  chamber  widths  of  42  and 
30mm,  and  less  than  3%  deviation  in  Nu,  for  all  other  chamber  widths.  These  maximum 


24 


deviations  corresponded  to  a  temperature  dilTerence  between  components  having  maxi¬ 
mum  and  minimum  average  temperatures  of  l.l^C  for  the  vertical  component  orien¬ 
tation,  and  6‘’C  for  the  horizontal  component  orientation. 

Aytar  (1991)  found  that  array  averaged  Nussclt  (Nu,)  and  flux  based  Rayleigh 
(Ra^)  numbers  plotted  as  straight  lines  of  nearly  constant  slope,  independent  of  chamber 
width,  provided  that  numbers  corresponding  to  0.1 15W  were  omitted.  He  reported 
a  range  of  slopes  between  0.296  and  0.340.  Using  a  similar  analysis,  Nu,  was  plotted  as 
a  function  of  the  temperature  based  Rayleigh  number  (Ra,)  and  chamber  width  for  both 
the  horizontal  component  orientation  (1  x  10‘<  Ra,<  2  x  lO’)  and  the  vertical  compo¬ 
nent  orientation  (3x  10’<3x  10*)  .  These  data  are  presented  in  Figure  9.  From 
Figure  9,  the  combined  data  set  forms  straight  lines  of  constant  slope,  independent  of 
chamber  width.  The  effect  of  chamber  width  on  the  Nusselt  number  is  more  pronounced 
at  18,  11  and  7mm,  whereas  at  42  and  30mm,  the  effect  is  almost  negligable. 
a.  Effect  of  Rayleigh  Number 

In  a  study  conducted  by  Aytar  (1991)  using  a  horizontal  component  orien¬ 
tation,  Rayleigh  number  dependence  was  determined  using  a  flux  based  Rayleigh  num¬ 
ber.  Prior  to  the  correlation  of  any  data,  Nusselt  number  was  plotted  as  a  function  of 
both  the  temperature  and  flux  based  Rayleigh  numbers,  using  the  data  sets  from  both 
the  horizontal  and  vertical  orientations.  Single  straight  line  behavior  was  observed  for 
both  data  sets  when  the  temperature  based  Rayleigh  number  was  used  (Figure  9). 
Based  on  this  observation,  it  was  determined  that  the  temperature  based  Rayleigh 
number  was  more  appropriate  than  the  flux  based  Rayleigh  number.  The  relationship 
between  the  Nusselt  number  (Nu,)  and  the  temperature  based  Rayleigh  number  {Ra,) 
was  assumed  to  be  of  the  simple  form: 

Nu,  =  Ci7?af' 


where  c,  and  are  constants. 

Because  of  the  similar  behavior  of  components  within  a  row,  row  averaged 
as  well  as  array  averaged  values  were  determined.  These  values  were  plotted  for  each 
spacing  and  dielectric  liquid  using  a  curve  fit  software  program  entitled 
"TABLECURVE"  which  fits  the  data  using  a  least  squares  regression  algorithm.  Curve 
fit  equations  for  array  averaged  values  of  Nu,  and  Ra,  for  each  width  are  presented  in 
Figure  10  through  Figure  14.  Curve  fit  equations  for  bottom,  middle  and  top  row  av¬ 
eraged  values  of  Nu,  and  Ra,  for  each  chamber  width  are  presented  in  Figure  15through 
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averaged  values,  since  presenting  correlations  for  all  of  the  rows,  as  well  as  the  array, 
would  be  lengthy  and  cumbersome.  However,  all  the  pertinent  data  is  contained  within 
to  produce  a  row  averaged  correlation,  should  the  need  exist. 

The  range  of  values  for  the  exponent  b^  from  the  array  averaged  data  for  the 
vertical  orientation  was:  0.352  <  b^  <  0.409,  with  an  average  value  of  0.381.  A  value  of 
0.35  was  chosen  to  represent  the  effect  of  Ra,  for  the  entire  array  at  all  chamber  widths. 
Aytar  (1991)  reported  using  a  value  of  =  0.30. 
b.  Effect  of  Chamber  Width 

The  effect  of  chamber  width  was  also  assumed  to  be  of  the  form: 


Nu,  =C2X^^ 


where  X  (~jy')  is  the  non-dimensional  chamber  width  and  W  is  the  maximum  chamber 
width  attainable  (42mm).  The  constant  Cj  reflects  the  Rayleigh  number  dependence. 
Using  the  curve  fit  equations  developed  for  each  chamber  width,  and  an  arbitrary  value 
of  Ra,  =  2  X  10’,  the  Nusselt  number  was  plotted  as  a  function  of  the  non-dimensional 
chamber  width  (see  Figure  30),  resulting  in  the  following: 

Nui  =  25.778X°-^^ 

and  therefore  the  following  general  correllation  for  FC-75: 

Nu,  =0.071 


which  is  valid  over  the  ranges; 

1  X  10*  <  <  3  X  10* 

23.04  <Pr<  30.73 

This  correlation  is  valid  for  low  values  oC  Ra,  (i.e.l  x  10^)  but  not  for  power 
levels  less  than  0.2W,  since  those  points  did  not  follow  the  trend  of  the  data  in 
Figure  9,  and  were  omitted  from  the  analysis.  This  correlation  is  accurate  to  within  6% 
of  the  curve  fit  equations  developed  by  TABLECURVE  from  the  original  array  averaged 
data. 

2.  FC-43 

The  maximum  deviation  in  the  flux  based  Rayleigh  number  (Ra^)  and  the 
Nusselt  number  (Nu,)  for  the  vertical  orientation  was  13.9%  and  16.9%  respectively. 
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Both  occurcd  at  a  spacing  of  11mm  and  a  power  setting  of  2.85W.  The  average  devi¬ 
ation  in  Ra^  and  Nu,  for  the  vertical  orientation  was  7.2%  and  12.2%  respectively. 
Aytar  (1991)  reported  an  approximate  deviation  of  12%  for  both  i?ay-and  Nui  at  a  7mm 
spacing.  These  ma.ximum  deviations  corresponded  to  a  temperature  difference  between 
the  components  having  the  maximum  and  minimum  average  temperatures,  of  7.4°C  for 
the  vertical  orientation,  and  8°C  for  the  horizontal  orientation. 

Similar  to  the  analysis  of  FC-75,  Nu,  was  plotted  as  a  function  of  Ra,  and 
chamber  width  for  both  the  horizontal  component  orientation  (3  x  10’  <  Ra,  <  5  x  10’) 
and  the  vertical  component  orientation  (8  x  10’</?a, <2x  10’).  These  arc  shown  in 
Figure  31.  Again,  data  points  corresponding  to  0.1 15W  were  omitted.  As  observed  for 
FC-75,  the  data  at  each  chamber  width  follows  a  straight  line,  with  nearly  constant 
slope,  independent  of  chamber  width.  Again,  the  effect  of  chamber  width  is  almost 
negligable  between  the  42  and  30nun  widths.  This  was  also  observed  by  Aytar  (1991). 
a.  Effect  of  Rayleigh  Number 

Similar  to  the  analysis  for  FC-75,  the  relationship  between  the  Nusselt 
number  (Nu,)  and  the  temperature  based  Rayleigh  number  {Ra,)  was  assumed  to  be  of 
the  form; 


Nu,  =  c,Raf’ 

where  Z),  and  c,  are  constants.  Again,  because  of  the  similar  behavior  of  components  in 
the  same  row,  row  averaged  as  well  as  array  averaged  values  of  Ra,  and  Nu,  were  de¬ 
termined.  Curve  fit  equations  for  array  averaged  values  of  Nu,  and  Ra,  are  presented  in 
Figure  32  through  Figure  36.  Curve  fit  equations  for  bottom,  middle  and  top  row  av¬ 
eraged  values  of  Nu,  and  Ra,  are  presented  in  Figure  37  through  Figure  51.  For  the 
same  reasons  given  previously,  only  array  averaged  values  were  used  after  this  point. 

The  range  of  values  for  the  exponent  from  the  array  averaged  data  for  the 
vertical  orientation  was;  0.354  <m<  0.387  ,  with  an  average  value  of  0.371.  A  value  of 
fc,  =  0.35  was  chosen  to  correspond  with  the  value  selected  for  FC-75. 
b.  Effect  of  Chamber  Width 

The  relationship  between  the  non-dimensional  chamber  width  and  Nusselt 
number  was  again  assumed  to  be  of  the  form; 

Nu,  =  CjX*' 
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Using  the  curve  fit  equations  developed  for  each  chamber  width,  and  an  arbitrary  value 
of  Ra,  =  2  X  10’,  the  Nusselt  number  was  plotted  as  a  function  of  the  non-dimensional 
chamber  width,  which  resulted  in  the  following  relation  (see  Figure  30): 

Nuj  =  0.0S03y?a?-^ 


which  is  valid  over  the  ranges: 

3  X  10^  </?«,<  1  X  10* 

57.3  <Pr<  127.5 

Similar  to  the  correlation  developed  for  FC-75,  this  correlation  is  valid  for  low  values 
of  Ra„  but  not  for  low'  power  levels,  since  those  data  points  w'ere  ommitted  from  the 
analysis.  This  correlation  is  accurate  to  within  3%  of  the  curve  fit  equations  derived 
from  the  array  averaged  data  using  TABLECL’RVE. 

3.  Effect  of  Component  Orientation 

The  effect  of  component  orientaion  on  the  heat  transfer  characteristics  of  the 
array  was  found  to  be  inherently  accounted  for  within  the  Grashof  number: 

Gr  = - 5 - 

V 

provided  that  the  dimension  of  the  component  in  the  direction  of  the  gravitational  vector 
(g)  was  used.  This  result  was  based  on  a  straight  line  fit  of  the  data  from  both  orien¬ 
tations,  when  the  length  scale  used  in  the  temperature  based  Rayleigh  number  was  cho¬ 
sen  in  the  direction  ofg. 

4.  Effect  of  Prandtl  Number 

As  in  the  other  cases,  the  relationship  between  the  Nusselt  number  and  the 
Prandtl  number  was  assumed  to  be: 

Nu,  = 

where  Cj  is  a  constant  which  reflects  the  dependence  on  the  Rayleigh  number  and  the 
non-dimensional  chamber  width.  In  this  analysis,  an  average  value  of  —  0.25,  the 
non-dimensional  chamber  w’idth  exponent  has  been  used.  Because  the  Prandtl  number 
for  each  dielectric  varied  significantly  over  the  range  of  chamber  widths  and  power  levels 
investigated,  an  average  value  of  Pr  was  used  for  each  dielectric  liquid.  Assuming  an 
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for  each  dielectric  varied  significantly  over  the  range  of  chamber  widths  and  power  levels 
investigated,  an  average  value  of  Pr  was  used  for  each  dielectric  liquid.  Assuming  an 
arbitrary'  value  of  Ra,  =  2x10’  ,  the  Nusselt  number  was  calculated  for  each  non- 
dimensional  chamber  width  and  average  Prandtl  number  using  the  general  correlations 
found  for  each  dielectric  liquid.  The  Nusselt  number  was  then  plotted  as  a  function  of 
average  Prandtl  number  for  the  various  non-dimensional  chamber  widths  (Figure  52). 
From  the  5  curves,  average  values  of  aandbi  were  determined.  This  led  to  the  overall 
correlation: 


which  is  valid  over  the  ranges: 

3  X  10^  <  <  3  X  10® 

23.04  <  Pr<  127.5 

This  correlation  is  accurate  to  within  1 1%  of  the  curve  fit  equations  developed  from  the 
original  array  averaged  data  for  FC-75  and  FC-43.  This  accuracy  has  been  determined 
using  the  full  range  of  Prandtl  numbers  for  both  dielectrics,  rather  than  simply  the  av¬ 
erage  values  used  to  determine  the  correllation. 

It  is  important  to  note  that  although  the  Prandtl  number  dependence  in  the 
overall  correlation  is  an  order  of  magnitude  less  than  either  the  Rayleigh  number  or 
chamber  width  effects,  only  2  data  points  for  the  Prandtl  number  were  used  to  determine 
this  dependence.  Additional  experimental  analysis  may  prove  that  this  dependence  is 
something  other  than  what  was  presented  here. 

An  interestmg  result  is  obtained  if  bi,  the  exponent  of  Ra„  is  chosen  as  0.33  or 
1/3.  The  motivation  for  this  choice  is  that  the  heat  transfer  coefficient  becomes  inde¬ 
pendent  of  the  components  length  scale,  as  seen  for  turbulent  flows.  If  the  analysis  used 
to  arrive  at  the  overall  correlation  presented  above  is  repeated  with  b^  =  0.33,  the  accu¬ 
racy  of  the  the  resulting  correlation  is  essentially  the  same  as  the  initial  correlation  de¬ 
rived  using  =  0.35.  The  correlation  which  resulted  from  this  analysis  was:,  ' 

Nu,  =0.1005/?a°®®X°-^®Pr°-°^® 
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5.  Uncertainty  Analysis 

Using  the  method  of  Kline  and  McClintock  (1956),  a  zeroth  order  uncertainty 
analysis  was  carried  our  for  all  chamber  widths  and  power  levels,  to  determine  the  un¬ 
certainty  in  the  Nusselt  number  (Nu,)  and  temperature  based  Rayleigh  number  (Ra,). 
The  calculations  were  done  on  a  component  basis.  The  expressions  used  in  this  analysis 
are  presented  in  Appendix  A.  Any  uncertainty  in  calculating  the  thermophysical  prop¬ 
erties  of  the  dielectrics  has  been  neglected.  Upon  reviewing  the  data,  the  following 
obersvations  were  made: 

1.  The  maximum  uncertainty  in  both  Nu,  and  Ra,  was  2.5%  for  FC-75  and  1.9%  for 
FC-43. 

2.  The  maximum  uncertainty  in  Ra,  and  Nu,  occured  at  the  lowest  power  level, 
0.115W.  This  held  true  for  both  dielectric  liquids.  This  result  was  expected  since 
the  ma.ximum  uncertainty  in  the  temperature  measurements  occured  at  this  power 
level. 

3.  The  uncertainty  in  Nu,  and  Ra,  decreased  slightly  as  the  chamber  width  decreased 
or  the  power  level  increased. 
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Figure  7.  Average  Array  Temperature  as  a  function  of  Net  Poner  and  Chamber 
Width  for  FC-75. 
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Figure  8.  Average  Array  Temperature  as  a  function  of  Net  Power  and  Chamber 
Widtii  for  FC-43. 
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Figure  9.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  all  Chamber 
Widths. 
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Figure  10.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  42mm 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  11.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  30mm 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  12.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  18mm 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  13.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  llnim 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  14.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  7ram 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  15.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  42mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  16.  Nu,  vs  Ra,  for  a  Vertical  Componetit  Orientation,  FC-75  and  a  42nitn 
Chamber  Widtli  using  Middle  Row  Averaged  Values. 
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Figure  17.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  42mm 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  18.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  30mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  19.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  30nim 
Chamber  Width  using  Middle  Row  Averaged  V'alues. 
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Figure  20.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  30mm 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  21.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  18nim 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  22.  NUj  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  18mm 
Chamber  Width  using  Middle  Ro>v  Averaged  Values. 
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Figure  23.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  18mni 
Chamber  Width  using  Top  Row  Averaged  Vaiues. 
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Figure  24.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  11mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  25.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  llmin 
Chamber  Width  using  Middle  Row  Averaged  Values. 
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Figure  26.  NUi  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  an  Ilnim 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  27.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  7inm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  28.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  7nim 
Chamber  Width  using  Middle  Row  Averaged  Values. 
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Figure  29.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-75  and  a  7mm 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  30.  Nu,  vs  X  for  a  Vertical  Component  Orientation  and  both  Dielectric 
Liquids. 
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Figure  31.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  all  Cham¬ 
ber  Widths. 
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Figure  32.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  42nim 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  33.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  30mni 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  34.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  18mm 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  35.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  llnini 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  36.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  7mm 
Chamber  Width  using  Array  Averaged  Values. 
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Figure  37.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  42mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  38,  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  42mm 
Chamber  Width  using  Middle  Row  Averaged  Values. 
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Figure  39.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  42mm 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  40.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  30mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  41.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  30mm 
Chamber  Width  using  Middle  Row  Averaged  Values. 
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Figure  42.  Nu,  vs  Ha,  for  a  Vertical  Component  Orientation,  FC-43  and  a  30ram 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  43.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  18mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  44.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  18mm 
Chamber  Width  using  Middle  Roav  Averaged  Values. 
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Figure  45.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  1 8mm 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  46.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  11mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 


70 


lOOOOO  1 000000  le+0?  le+08  le+09 

Rf=^YLEIGH  NUMBER 


Figure  47.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  11mm 
Chamber  Width  using  Middle  Row  Averaged  Values. 
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Figure  48.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  an  Ilmni 
Chamber  Width  using  Top  Row  Averaged  Values. 
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Figure  49.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  7mm 
Chamber  Width  using  Bottom  Row  Averaged  Values. 
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Figure  50.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  7mm 
Chamber  Width  using  Middle  Row  Averaged  Values. 
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Figure  51.  Nu,  vs  Ra,  for  a  Vertical  Component  Orientation,  FC-43  and  a  7mm 
Chamber  Width  using  Top  Ron  Averaged  Values. 
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Figure  52.  Nui  vs  and  Non-Dimensional  Chamber  Width  for  a  Vertical  Com¬ 
ponent  Orientation. 
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VI.  CONCLUSIONS 


The  natural  convection  heat  transfer  characteristics  of  a  3x3  array  of  vertically  ori¬ 
ented  simulated  electronic  components  were  studied.  The  components  were  mounted  to 
a  vertical  plate  affixed  within  a  rectangular  enclosure  filled  with  a  dielectric  liquid  (FC-75 
or  FC-43).  Power  was  supplied  to  the  components  using  Inconel  foil  heaters,  and  varied 
from  0.1 15W  to  2.90W  per  component.  The  chamber  width,  measured  from  the  surface 
of  the  circuit  board  to  the  opposite  enclosure  wall,  varied  from  42mm  to  7mm.  This  was 
accomplished  using  plexiglass  spacers.  Data  from  this  study  was  correlated  with  the 
data  obtained  from  a  horizontal  component  orientation,  which  was  investigated  by 
Aytar  (1991).  The  following  conclusions  are  made; 

1.  A  general  emperical  correlation  which  predicts  the  Nusselt  number  as  a  function 
of  temperature  based  Rayleigh  number,  non-dimensional  chamber  width,  Prandtl 
number  and  component  orientation  was  determined.  This  correlation  was  based 
on  array  averaged  values  for  the  Nusselt  and  Rayleigh  numbers. 

Nui  =  0.0717/?a°-^ 

3  X  10^  <  <  3  X  10* 

23.04  <Pr<  127.5 

The  accuracy  of  this  correlation  was  within  11%  of  the  original  array  averaged 
data. 

2.  The  average  heat  transfer  coefficient  is  independent  of  the  component  length  scale 
if  bi,  the  Rayleigh  number  exponent,  is  chosen  as  0.33.  A  correlation  was  deter¬ 
mined  that  provided  results  with  nearly  the  same  accuracy  as  the  correlation  given 
above. 

3.  The  effect  of  component  orientation  was  found  to  be  inherently  accounted  for 
within  the  temperature  based  Rayleigh  number  {Ra),  provided  the  dimension  of  the 
component  in  the  direction  of  the  gravitational  vector  (g)  was  used  as  the  length 
scale  in  Ra, . 

4.  The  effect  of  chamber  width  on  the  Nusselt  number  was  almost  negligable  between 
the  42nim  and  30mm  chamber  widths.  However,  as  chamber  width  varied  from 
18mm  to  7mm,  its  effect  on  the  Nusselt  number  became  more  pronounced. 

5.  Heat  transfer  in  FC-75  was  largely  due  to  convection  resulting  from  buoyancy 
forces,  regardless  of  chamber  width. 

6.  Heat  transfer  in  FC-43  at  chamber  widths  of  11mm  or  greater  was  mostly  due  to 
molecular  diffusion. 
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7.  Heat  transfer  in  FC-43  at  a  chamber  width  of  7nim,  in  contrast  to  the  results  for 
the  other  chamber  widths,  was  due  mostly  to  convection  resulting  from  buoyancy 
forces. 

8.  The  maximum  average  component  temperature  for  the  horizontal  orientation  was 
7®C  higher  than  the  vertical  orientation,  when  FC-75  was  used.  For  FC-43,  the 
maximum  average  component  temperatures  were  nearly  identical  for  both  orien¬ 
tations. 

9.  At  a  maximum  heat  flux  of  approximately  17  Watts/cm’,  the  highest  average  com¬ 
ponent  temperature  for  the  vertical  orientation  was  more  than  20‘’C  below  the 
typical  maximum  allowable  component  temperature  of  85®C. 

10.  The  maximum  uncertainty  in  the  Nusselt  and  Rayleigh  numbers  was  2.5%  and 
occured  at  the  lowest  power  level  investigated  (0.1 15W  per  component). 
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VII.  RECOMMENDATIONS 


Based  on  the  findings  of  this  study,  the  following  recommendations  are  made: 

1.  Investigate  additional  dielectric  liquids  with  range  of  Prandtl  numbers  wide  enough 
to  determine  a  more  accurate  Prandtl  number  dependence. 

2.  Mount  thermocouples  to  the  substrate  in  order  to  account  for  heat  loss  from  the 
substrate  to  the  dielectric  liquid. 

3.  Correlate  the  horizontal  and  vertical  orientation  data  using  a  more  elaborate  form, 
in  order  to  improve  the  accuracy  of  the  Nusselt  number  prediction. 

4.  Investigate  a  variable  geometry,  non-uniform  component  orientation. 
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APPENDIX  A.  UNCERTAINTY  ANALYSIS 


Using  the  method  of  Kline  and  McClintock  (1956),  a  zeroth  order  uncertainty 
analysis  was  carried  out  to  determine  the  uncertainties  in  the  Nusselt  number  (Nu,)  and 
the  temperature  based  Rayleigh  number(i?a,),  for  each  of  the  chamber  widths  and  power 
levels  investigated.  The  uncertainties  in  Nu,  and  Ra,  were  calculated  on  a  component 
basis.  Any  uncertainty  in  calculating  the  thermophysical  properties  of  the  dielectric 
liquids  has  been  neglected.  The  expressions  used  in  this  analysis  are  presented  below. 
Numerical  values  have  been  omitted  to  maintain  generality, 

1 .  Uncertainty  in  the  power  supplied  to  the  heating  elements 

(volt  —  emf) 

POWER  =  emf- - j ^  =Aemfvoll,  R^) 


SPOWER  = 
where; 


d  POWER 

demf 


bemfy  + 


POWER 

dvolt 


bvoll 


b  POWER  (voli-lemf) 


bemf 

b  POWER 

bvoll 


Pp 

emf 

Rn 


d  POWER 

dRr. 


bR. 


bpower  (volt  —  emf 

=  -  emf- 


bR„ 


R 
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2.  Uncertainty  in  the  heat  loss  through  the  plexiglass  substrate 

Cte,  =  =/Ar,Rj 


^Qloss  = 


where: 


'OSS 


BAT 


bATY  + 


SQu 


OSS 


bR, 


•  bR, 


bemf  =  bvolt  =  +  0.001  V 
bRp  =  ±  0.05Q 


^Qloss 

bAT 


R. 
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I^I2£L  =  _AL 
SRc  Ri 

SAT=  10 
SR,=  10 

3.  Uncertainty  in  the  net  heat  dissipated  by  the  heaters 

=  POWER  -  Qi,,,  =APOWER,Qi,,,) 


SQne,  =  sJiSPOWERf  +  iSQi,J 


4.  Uncertainty  in  the  average  heat  transfer  coefficient 

Qnet 


h  = 


^  total' 
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J^Qnel,A  total, ^  P) 
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A,^,^i  =  2BH+2LH+B 

total 
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total 


dL 


dB 

total 

dH 

total 

dL 


=  2H+L 
=  2H+  B 
=  2H+  B 


SAT=  1 

SB  =  SH=SL=  \  x  10"^»i 


5.  Uncertainty  in  the  Nusselt  number 


Nu  = 


hL 
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where: 


axu  _  L 

dh  kf 

aNu _ h_ 

dL  ~  kf 


6.  Uncertainty  in  the  Temperature  based  Rayleigh  number 


Ra,  = 


va 


=AL,M) 


5  Raj 


+ 


where: 

dRa, 

8L  ~ 

cRa,  _  gfiL^ 
BAT  ~ 
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